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Abstract 

The (^-function was introduced by Affleck and Ludwig in the context of critical 
quantum systems with boundaries. In the framework of the thermodynamic Bethe 
ansatz (TBA) method for relativistic scattering theories, all attempts to write an 
exact integral equation for the off-critical version of this quantity have, up to now, 
been unsuccesful. We tackle this problem by using an n-particle cluster expansion, 
close in spirit to form-factor calculations of correlators on the plane. The leading 
contribution already disagrees with all previous proposals, but a study of this and 
subsequent terms allows us to deduce an exact infrared expansion for g, written 
purely in terms of TBA pseudoenergies. Although we only treat the thermally- 
perturbed Ising and the scaling Lee- Yang models in detail, we propose a general 
formula for g which should be valid for any model with entirely diagonal scattering. 
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1 Introduction 



The study of two-dimensional conformal field theories with boundaries 1_ and their 
integrable perturbations |21 El 15 is of interest both in condensed matter physics [S] 
and in string theory An important quantity emerging from the definition of the 
cylinder partition function for these field theories is the ^-function, the 'ground-state 
degeneracy' or 'boundary entropy', which for models critical in the bulk was introduced 
some years ago by Affleck and Ludwig 0|. While many interesting questions remain 
in these cases jSl |H] , in this paper we shall deal with the further issues which arise for 
off-critical, massive, systems. 

The 5- function for massive field theories can be defined as follows jlUl 1111 1121 113j . 
There are two possible Hamiltonian descriptions of the cylinder partition function. In 
the so-called L-channel representation the role of time is taken by L, the circumference 
of the circle: 

00 

Za, = TVH(.,,,e-^<r(^-^-«) = Y: e-^^""^(*^'^) . (1.1) 

n=0 

In this formula, H^^^^ propagates states in 'H(^a,(})^ the Hilbert space for an interval 

of length R with boundary conditions a and (3 imposed at the two ends, En^^^ G 
spec(-ff^*^^'') , and M is the mass of the lightest particle in the theory. In the R-channel 
representation the role of time is instead taken by R, the length of the cylinder: 



Z,^ = (a|e-«^"^^(^'^)|/?) = g(")(/)g^'^)(Oe-^^""^(^'^), {l = ML) (1.2) 

n=0 

where E"'^^ G spec(if"'''^) and 

Qt\l) = j^^.. (1.3) 

In equation (|1.2|) . the boundary states |a), \j3) and the eigenbasis {IV'n)} of the Hamil- 
tonian H™'^ have been used. These are defined on a circle of circumference L and 
propagate along the 'time' direction R. At large Z, the function ln^i''''(/) grows linearly: 

lngW(/)~-/„L, (1.4) 

where the constant fa contributes to the constant (boundary) part of the ground-state 
energy on the strip (see eq. (jA.Si) ). The standard (/-function is then defined as 

\nga{l)=\ng^^\l) + faL . (1.5) 

In theories with only massive excitations in the bulk, \nga{l) tends exponentially to 
zero at large /. 
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R 

(space) 

Figure n The L-channel decomposition; states \xn) live on the dotted hne segment 
along the cyhnder. 




R 

(time) 

FigureEl The R-channel decomposition; states \ipn) live on the dotted circle around 
the cylinder. 

The two decompositions are illustrated in figures ^ and l2j 

The equality of and results in the following important identity: 

oo oo 
n=0 n=0 

The purpose of this paper is to develop an exact expression for the ground-state function 
through the large- -R limit of (|1.6() . with boundary conditions (3 = a. As it 
stands, the fact that E™^{M, L) is negative makes the RHS of (|1.6|) diverge as i? ^ oo; 
however, rearranging gives 

21ng(°)(/) = REf^{M,L)-LEf'''{M,R) (1.7) 

+ In ^1 + ^ ^-L{ET"'{M,R)-Ef"^{M,R))\^ ^ ^^^_K(Eci.c_sci.c) ^ ^ 

We shall restrict our attention to massive theories with non-degenerate ground state 
on the plane. For these models the non-zero mass gap gives the final term the leading 
behaviour 

0(g-ij(£;f-(M,L)-i?S'-(M,L))) _ 0{e-^M) (1.8) 
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in the domain R ^ L ^ 0. In this same domain, "^[M, R) tends to its hmiting form 
as 

Ef'^{M, R) = SM^R + 2U + O(e-^^0 (1.9) 

where £ and fa are the extensive bulk and boundary free energies, as in (TO) . These 
constraints are crucial for the validity of the perturbative treatment to be introduced 
in the following sections: the higher corrections have a clear dependence on R and do 
not contribute to the (jf-function. Discarding these exponentially-suppressed terms and 
using the definition (jl.Sj) . we finally obtain 

21n<7i°)(/) ~ R[Et%M, L) - SM^l) + e-^(i^r-{M.R)-<"^(M,i?))^ . 

(1.10) 

Having taken R to be large, the cluster expansion involves letting L tend to infinity 
as well, so that an expansion of the RHS of (|1.10j) can be developed in terms of one-, 
two- and so on particle contributions, which themselves can be consistently estimated 
using the Bethe-ansatz approximated levels (|A.6|) , HA.7|) . Note that this differs from the 
strategy adopted in ^01 ) where a saddle-point evaluation of the dominant contributions 
at finite L was made instead. 

The rest of this paper is organised as follows. In section |2l the cluster method is 
exemplified by studying the free fermion theory associated to the thermally-perturbed 
Ising model. The resulting integral expression for \D.ga{l) turns out to be in full agree- 
ment with previous results of ^1 . In section |51 two previous proposals ^1 E] 
for \n.ga{l) are described and in section |3 the scattering data for the scaling Lee- Yang 
model, our working interactive example, are summarised. The ultraviolet result ob- 
tained from the conformal perturbation theory and the boundary truncated conformal 
space approximation (BTCSA) jl2| ll,S| IH] is compared with infrared numerical results 
from the Bethe Ansatz, and the equivalence between the two functions is confirmed by 
a large overlap at intermediate scales. 

This agreement motivates the search for an exact analytic expression. This is the 
main objective of sections 15. H 15.21 and 15.31 where the large strip-width (i? — > oo) limit 
is explicitly taken and sums over the quantum numbers are transformed into integrals 
in rapidity variables. This analysis leads to the final exact expansion for \n.gci{l) given 
in equation (|5.24|) . This and its generalization (|5.30() constitute the main results of the 
paper. In section l5^ we also briefly comment on the similarity between our results and 
one recently obtained by Woynarovich in ^Sj. Section |^ contains our conclusions. Fi- 
nally, in Appendix^we summarise the main equations used to develop our programme: 
the thermodynamic Bethe ansatz ^H] and the Bethe quantisation conditions. In Ap- 
pendix^ the reflection factors for the Ising model are recalled and explicit expressions 
for the boundary entropy for free-free and fixed-fixed conditions are presented. 
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2 A simple example: the Ising model 

We start with the study of the free Majorana fermion theory corresponding to the 
thermahy-perturbed Ising model on a strip. The partition function is (cf. * : 

Z„« = e-i^-^CM.fl) + e-^coshe,^ ^ ^2.1) 

i>o 

or 

lnZ„^ = -L^^*"P(M,i?) + J^ln(l + e-'™'^'^^0 

j>0 

1 °° 1 
= -LEf"'{M,R) + - ln(l + e-'™^^^0-^ln(l + e-'). (2.2) 



2 

j=-oo 



Due to the singular behaviour for the Ising model of the bulk and linear terms £ and fa 
defined in (|A.2|) and HA.5|) below, it is convenient, exceptionally for this case, to work 
with subtracted energies tending exponentially to zero at large scales: 

Ef''P(M,R) ~0 (2.3) 

and consistently to set 

E^'^{M,L) = - / — Mcoshe ln(l + 6-'™*^^^) (2.4) 

(cf. ()A.2|) '). Starting from the quantization condition 

r sinh % - i In Ra{0j) = vrj , (r = MR) (2.5) 

with integer j and Ra{Oj) as defined in HB.1|) . writing equation (|2.5jl with J — > J + 1 and 
subtracting (|2.5|1 from the result, we find in the large R limit 

^ {rcosh{dj)-i^\nR^{6j)^ +0{{Mjf) = 1. (2.6) 

Substituting this into (|2.2|) . 

lnZ«„~^ / (i0 f - cosh(0) + 0^(0)- 5(0)) ln(l + e-'™^^^), (2.7) 

where (pai^) is given in HB.2|) . In the latter equation we recognize a part corresponding 
to REq"^{M,L), and, considering also (|2.3|1 . we arrive at the exact result 



21n5„(/) = hm {In Zaa{L,R)+RE^''''{M,L)) = ]- f dO {^^{0) - 5{e))Hl+e 



■I cosh 



'Notice that the zero momentum {6j — -i-^ j — 0) particle state is forbidden. 



4 



This coincides with the result found in |1U[ lllj using a different technique. The match 
confirms the correctness of our method, at least in this case, and motivates its study in 
more complicated models. 

In figures Ol and \^ the integration of ()2.8|) for free-free and fixed-fixed boundary 
conditions, corresponding to k = 1 and A; = — oo in ()B.1|) and (|B.2|) ^. is compared with 
numerical results obtained by estimating the large- i? partition function ()2.1|) using the 
Bethe ansatz quantized energy levels (|2.5|) directly, and then extracting the boundary 
entropy through the relation 



21n<7a(0 ~ {\nZ^^{L,R)+RE^"'{M,L)) 



(2.9) 




Figure 121 2 In gfixed vs. ln(Z) for Ising with 
fixed boundary conditions. From the bottom, 
the lines represent 5, 100, 500, and 5000 
particle contributions. The maximum quantum 
number is 80 and r = 10. The top line is the 
exact result. 



Figure^ 2 In gtree vs. ln(Z) for Ising with free 
boundary conditions. From the bottom, the 
lines represent 5, 100, 500, and 5000 particle 
contributions. The maximum quantum number 
is 80 and r = 10. The top line is the exact 
result. 



For interacting models a compact expression such as 1)2. 1|) is not available, and one 
is forced to build the partition function using the LHS of (|1.6j) directly. In figures El 
and we test this more general way to estimate a g(- function. A similar idea was first 
applied to the scaling Lee- Yang model in |12|ll3j: however in that case the energy levels 
were estimated using the BTCSA method ^^l; rather than the Bethe ansatz. 

^By studying the monodromies of the integral 112. 8|l . we have also found more explicit expressions for 
ln5flxed(0 ^iid ln(?free(Oi these are given in Appendix lEl 
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-1.5 -1 -0.5 0,5 1 1.5 2 

Figure ini 2 In gflxed vs. ln(/) for Ising with 
fixed boundary conditions. From the bottom, 
the lines represent cluster contributions of 
1, 2, . . . , 8 particles. The maximum quantum 
number 80 and r = 10. The top line is the 
exact result. 



-2.5 -2 -1,5 -1 -0.5 0.5 1 1.5 2 

Figure El 2 In g^ree vs. ln(Z) for Ising with free 
boundary conditions. From the bottom, the 
lines represent cluster contributions of 
1, 2, . . . , 8 particles. The maximum quantum 
number 80 and r = 10. The top line is the 
exact result. 



We see that it is hard to get a good estimate of the ultraviolet value of In g from this 
form of the cluster expansion. In section 0] we shall solve this numerical problem for the 
case of the scaling Lee- Yang model by matching this numerical Bethe ansatz calculation 
with the ultraviolet perturbed CFT results of ^Hl, while in section El we shall develop a 
more analytical treatment. 



3 Earlier proposals for g 

Consider a 1+1 dimensional integrable field theory with entirely diagonal scattering and 
particle species. According to the proposal of TUI the boundary entropy should be 
given by an expression of the form 

1 ^ r 

Ing^il) = -,y^ de Qa{e) ln(l + e'^-^W) , (3.1) 
^ a=i 

where the function £a{S) is the solution of the periodic-boundary-conditions TBA ()A.H) . 
and 

Qa{0) = (0(^^(0) - 2<l)aa{2e) - 5(0)) (3.2) 

with 

<t^^a\0) = -lAin4«)(0) , <^,,(0) = -l-±l^Sam ■ (3.3) 
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(Note, the normalisations of (l)a (0) and 4>ab{(^) differ from those in by factors 

of vr and 2it respectively. This change is merely to simplify some later formulae.) 

However, the detailed analysis of showed that, for non-zero values of the lightest 
bulk mass M, the resulting /-dependence was incorrect, both in the total change in 
Qai^) between UV and IR, and in the behaviour of the small-/ series expansion. On the 
other hand, the predictions of 1)3. 1() and ()3.2|) for dependence of ga{l) on the boundary 
parameters at fixed /, and also for the ratios of (^-functions ga{l) / gp{l) , were in very 
good agreement with conformal perturbation theory and the BTCSA. This suggested 
that the formulae should be modified by some boundary condition independent extra 
terms, but provided little clue as to what those extra terms should be. 

Subsequently, it was proposed in that (|3.2j) should be replaced by 

Qa{0) = [(t>^^HO) - 24,aa{2e) - <Paa{0)) . (3.4) 

However, using results tabulated in [IB] it can be checked that this modification does 
not cure the problems arising in the bulk-massive case. 



4 The scaling Lee- Yang model 

The spectrum in the bulk scaling Lee- Yang theory consists of a single particle species, 
with two-particle bulk scattering amplitude ^1 IT^ 

S{e) = -{1){2) , (x)= .6 . (4.1) 

When a boundary is present, two different types of boundary conditions arise, which 
were labelled 1 and <I>(/i) in ^12 . The corresponding reflection factors are 

Rm(^) = Rb{e) , Riie) = Roio) , (4.2) 

where 

h \hc\ sm{{b + .5)7r/5)M^/^ , he = -0.6852899 . . . (4.3) 

is the coupling of the boundary field and 

Rm = il) (i) (I)"' {S{e + in'-^)S{9 - ivr^))-^ . (4.4) 

We first use the Bethe ansatz equation together with (|1.7() to obtain the (7-function up to 
six-particle contributions. The results are shown in figure |7| for the boundary condition 
1, and are compared with the ultraviolet expansion obtained from (boundary) conformal 
perturbation theory and the BTCSA jl3j : 



2 In ,1 (0 = i In + 2 ^ ^ + 2 E d„ (^) 0(/-) , (4.5) 
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where = i(^/3 - 1)M, di w -0.25312, da ^ 0.0775, « -0.0360, (i4 0.0195, and 
K PS 2.6429. 




FigureH 2 In vs. ln(/) with (1, 1) boundary conditions. The dotted hne is the 4**^ 
order result of ^31 and the sohd lines are BA results with r = 8. The solid lines, from the 
bottom, represent the 1, 2, • • • , 6 particle contributions respectively. The maximum 
quantum number used was 80 (57 for the 6-particle contribution). 

As can be seen from the figure, the results from the Bethe ansatz and from perturbed 
conformal field theory overlap over a significant range of scales. This supports our 
hypothesis that the two approaches are describing the same function g^{l), expanded 
about either the IR or the UV. 

5 Infrared expansion for the Lee- Yang model 

The purpose of this section is to develop an analytic technique to check the earlier 
proposals described in section|21and at the same time to give hints about the appropriate 
modifications. The idea, successfully applied above to the Ising model, is to start from 
the Bethe ansatz and to set up a cluster expansion by transforming the sums into 
integrals as ii ^ oo. The method is quite powerful, and already at first order it confirms 
the questions raised in [T^ about the proposals described in section |31 To simplify the 
discussion we shall only treat the (1,1) boundary conditions directly. However, this 
restriction is of no real significance since the results in 0J 1131 119j show that 

<i)«=^(-^)^rW' m=e^(^ (5.1) 

where e{6) is the solution of the ground state TBA equation with periodic boundary 
conditions. 
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5.1 The one particle contribution 

We start from the large- i? equation (|1.1U|1 . truncated at the one-particle level: 

2 In 51 r(^E^'''{M,L)-£M^L^ +ln(l+ ^ g-Zcoshei(ni) ^ _ ^ ^5 2) 

ni>0 

where the one-particle Bethe ansatz essentially coincides with that for a free particle, 

r sinhOi — iln Rj{9i) = nni , (5-3) 

with integer rii. Performing the continuous limit as in section |21 we find 

^ oo 1 /■ 

Pi = V e-'™^^^^i = L( y e^'™*^^^! -e-') — >- d0 - (5(0)) e"'"""^^ 

— ^ 2 \ ^ — ^ / 2 /lu 

ni>0 ni=-oo 

(5.4) 

where the Jacobian for the change of variable ni ^ 9i = 9 is 

j(i)(0) = -cosh + 0^(0). (5.5) 
vr 

The cosh term cancels the leading part of the term linear in R on the RHS of H5.2p . 
leaving the first contribution to In as 

21n<7i = ^ / d9{cl>j{9)-S{9))e-''°'''' + ... . (5.6) 

Comparing this result with the proposals of sectional we conclude that both are incor- 
rect: in particular no (j){29) or terms are involved in the leading large / asymptotic. 

Next, we want to use this result to gain a hint as to how the earlier proposals should 
be modified. However, the task to totally or even partially re-sum the cluster expansion 
directly is, in principle, very hard. Our work is driven by the extra assumption that 
the final result should depend, just like the earlier 'partially correct' proposals (|3.2j) 
and ()3.4() . on the bare single-particle energies only through the TBA pseudoenergies 
e{9). As will be reported in more detail below, the consistency of this assumption was 
checked carefully up to four particles and confirmed, by a more superficial inspection, 
to all orders. 

The attempt to find an exact expression for In^i, therefore, naturally starts from 

In5i = [ln5i]g^ + ... (5.7) 
where we have defined the 'dressed' version of (|5.6|) to be 

2[\ng^i^^ = \ [ (01(0)- 5(0)) ln(l+e-W). (5.8) 
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Figure ISl gives some initial numerical support for the conjecture. 




FigurelHl 2 In vs. In(^). The top dotted line is the 'exact' result obtained by combining 
the CFT and BA results from figure |3 The solid lines are obtained from H5.6|l (bottom 
line) and (|5.7|l (middle line). 

Notice now that [In (7]]^]^^ also contains n(> l)-particle contributions. To see this at 
second order, we expand 



ln(l + e 



as e 



-2e{e) 



+ ... 



(5.9) 



and use the exponential of equation HA.1|1 expanded in terms of the bare particle energy 
e{9) = M cosh (9 



g-e(e) _ g-Zcoshfr j ^ 



+ 



(5.10) 



to see that 

2[ln5l]g^ = ^ / dO {<t^^{e)-5{e))e- 

^ JR 

-] [ d0 0i(e)e-2'™^h^ + i / d^id^a 01(^1)0(^1 -^2)6-'^°^*^'^-'™^^'^ +... (5.11) 
4 Jm ^ Jr2 



The aim of the analysis in the next sections is to justify the replacement of g-'coshe 
by ln(l + e~^^^^) and is also to find some hints as to the origin and form of the further 
correction terms in 1)5. 7(1 . 
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5.2 Two and three particle contributions 

We start again from (|1.1UI) . this time keeping both one and two particle contributions: 



21n5i = r(^E^''%M,L) -SM'^L^ 

+ In (l + ^ g-zcoshei + ^ ^ cosh ei-z cosh 62 + . . ^ (5.12) 

ni>0 ni>0n2>ni 

where the two-particle-state momenta (^1,^2) are related to their quantum numbers 
(ni,n2) via the Bethe ansatz equations 

-sinhOi- -InRjiOi)- ^lnS{ei-e2)S{9i + e2) = m; 
IT vr zvr 

-sinh^2--lni?ii(^2)-7^1n5(02-ei)5(6'2 + 0i) = n2 . (5.13) 
TT vr zvr 

The new piece in ()5.12() can be written as 

_ ^ ^ g-; cosh 61 -« cosh 6*2 = i ^ ^ g-« cosh Si-Z cosh 6*2 

ni>0n2>"l ^ ni=— 00 n2=— 00 

00 



_ ^ g-Zcosh6»i-/ _ i ^ g-2Zcosh6»i _^ £g-2; _ (5.14) 
ni=— 00 ni=— 00 ^ 

As i? ^ 00 the continuous limit can be taken as: 

y y e-^""'^^'-^""'^^^ — >- [ deid02jf^(^i,^2)e-''°'^^i-'™'''^2 

8 /ll}2 

ni>0n2>ni 



/ d^jf (e)e-'™«^''-' - ^ / d^jf (0)e-2'-^t^^ + §e-2' (5.15) 

f2~) f2) f2~) 

The Jacobians Jr(ei,^2), ^2 (^) and J^^'{Q) can be calculated from the Bethe ansatz 
equations ()5.13|1 as before. Notice that the correct subtractions of the excluded contri- 
butions (those excluded by the statistics) are crucial to get the corresponding Jacobians 

f2) f2) 

J)^'{e) and J\ '{e): one has to take the derivatives only after the forbidden values of 

the quantum numbers n\ and are fixed. Although we have performed the calculation 

in full and checked the consistent cancellations of the r (strip size) dependent parts, for 

(2) (2) (2) 

brevity we shall concentrate on the subleading, r-independent, parts , ^2 ' is ■ 
2l\QxM) = </'l(^i)</'l(^2) + 20(0i- Wl(^2) + 20(0i-02)</'i(0i) 

= H{(^) + 2m, 3f\0) = H{(^) + 2<i^{2e). (5.16) 
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Expanding the logarithm in (|5.12|1 . we have at second order 

ln(l + Pi + P2 + . . .) = A + (^^2 - ^) + . . . . (5.17) 

The up-to second order 2lng-^ contains seven distinct contributions, the first five co- 
inciding with those written exphcitly on the RHS of (|5.11|) and corresponding to the 
up-to- two particle expansion of 2 [In gj]^^^- (This confirms the correctness of the dressed 
formulae H5.71 15. 8|) up to this point). There are also two genuinely new terms, and we 
find: 

21n5i = 2[ln5i]g) + i /" d^id^a ^ + W(^i " ^2)6"'^°^'^'^-'™^^'^ 

-If 0(20)6-2'™^^'^ + ... , (5.18) 

The final step is to iterate the dressing procedure, though in a modified form, by re- 
placing /coshS with e{9) and writing Ingi^ = [Ingj]^^^ + [ln(7][]^^ -|- . . . with 

2[ln<7l]Sp'^ = J / deid92ct>{ei+e2)m-02)e-'^''^-'^''^-l- [ d9cp{29)e-^'^'K (5.19) 

Again this dressing prescription can be justified retrospectively by testing at third and 
higher order. The third-order result turns out to support the assumption, and gives a 
genuinely new type of correction to llng^, independent of </'][(^): 



21n5i = 2[ln5i]i^^+2[ln5i]g^ + ^ /d^</>(20) 

"J Jr 



g 3( cosh t 



+ \ [ d0id02'ie3'/'(^l+^2)</'(^2-^3)'/'(^3-^l)e-''=°^^'^-'™^^'^-'^°^^'^ 

- / (i0id02</'(^i + 02)</'(^i -02)e~''°'^^''^-2'™'^^^ + ... . (5.20) 
5.3 The exact result 

To go further in the expansion becomes increasingly difficult, due to higher number of 
Jacobians and the huge number of terms contributing to a single Jacobian. However we 
managed to complete the analysis up to four particles and to perform a more superficial 
inspection at higher orders. The following results were deduced: at each order there is 
always a new contribution of the form 

Cn = - I dOi... d9n m + ^2)0(^2 " ^3) . • • 0(^n " ^1)6-'^°^^^^ . . . e"'™^^^" , (5.21) 

when the lower order terms C2, . . . , C„_i are corrected according to the rule 

g-/coshe^ ^ ^ /5_22) 

1 + e^(^) 
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An additional term contains a single integration over the function (f){29), as already 
seen in 1)5. 18() and H5.20() . To treat these terms is trickier, but after a few attempts we 
convinced ourselves that the correct procedure is always to replace I cosh 9 with £{9), 
and then to resum, as: 

1^ d9 (^-ie-2'-h. ^^29) ( ln(l + e'^^) - j^^) • (5-23) 

We finally arrive at 

21n5i(/) = ^ / (01 5(0)- 20(20)) ln(l + e-^W) 
°° 1 /" fl9 ii9 

+ E-/ T— 4^- --T-^ 0(^1+^(^2 -03)... 0(^n-0n+l), (5.24) 
n=l 

with 0n+i = 01 -f. It is not hard to check that this formula is consistent with the one-, 
two-, and three- particle results described above. Another simple (but nontrivial) check 
can be performed in the ultraviolet, as 

e=(e) ^ = ^ f d9cPj (9) = -2 , (5.25) 

^ Jm. 

[ d9m = -^, I d0i... d0„ 0(01 + 02)0(^2- ^3). •• 0(en - ^i) = ^4^, (5.26) 
and 1)5. 24(1 reduces to 

^^^^ ( 1 \ 1 11 

-ln(l + e--) + E^;^(TT7^)" = -Ml + e--)--ln(l + 



2n V 1 e-^o / ' 2 V 1 e^o 

n=l 

(5.27) 



2 V 2\/5 

which is the expected / = value for 2 In given in ()4.5p . 

The first term on the RHS of (|5.24)1 coincides with the proposal of 10,, repeated 
in equations 1)3. 1() and ()3.2|) above, while the remaining parts constitute a boundary- 
condition independent correction, which only comes nontrivially into play when the bulk 
theory is massive. As mentioned at the end of section |3J this was only to be expected 
given the results of |13j , but it is nevertheless satisfying that the already- verified portions 
of the earlier results have been recovered by this rather different route. 

In figure El the result from equation (|5.24|) is compared with the 'exact' result ob- 
tained in section |1] by combining UV perturbed CFT results with the IR cluster expan- 
sion from the Bethe ansatz. We see that keeping only the first three terms in the series 

* Notice the sequence of signs +,—,—,...,— in the arguments of the (^s in 15.241 . that the number 
of factors in the correction term is n and that according to Ij5.2i-i|l at n = 1 only (l>{26) survives. 
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already gives a very good agreement with the exact result: in the UV the agreement is 
within about 0.3%. 




FigiireEl 2 In gj^ vs. ln(/). The dotted line corresponds to the 'exact' contribution. The 
bottom solid line represents the RHS of H5.24|l with the sum truncated at the first term. 
The solid line just above the dotted line truncates the sum at the doubly-integrated term, 
and the line just below the dotted line is the total contribution up to the triple integral. 

We also used (|5.24() to make a numerical estimate of the coefficients of the ultraviolet 
expansion of lngj{l): 



lng-i{l) 0.3214826953191634 + 0.483692443734693x12 - 0.253117570 x 

„2 



0.037284 + 0.02042 x"^ 



0.0120 x^ + 0.0074 x^ + 



+ 0.0781176 X 
with X = (//k)12/5_ rpj^jg 

can be compared with the result of P^ : 

In gj{l) 0.3214826953191634 + 0.4836924437346968 x^ - 0.253117581 x 

+ 0.0775 x^ - 0.036 x^ + 0.0195 x"^ + . . . 



(5.28) 



(5.29) 



(Exact expressions for the first three coefficients in (|5.29j) were found in pSl ! here we only 
quote sufficiently-many digits to enable the numerical errors in ()5.28() to be assessed.) 

By inspection, it is straightforward to generalize (|5.24j) to more general theories 
with particle species and entirely diagonal scattering and reflection matrices. The 
proposal is 



1 ^ r 

2lngUl) = -Y, de{ct>^^\e)-5{e)-2ct>aam)Hl + e- 



■ea{e) 



+ 



oo TV 

E E 

n=l a\...an 



1 



de^ 



d0„ 



^ n Jjj" 1 + e^°i(^i) 



1 + e^a„{0n) 



Ui)) , (5.30) 
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where 6n+i = Oi, a„+i = oi, while (pabiO) and (pa are defined in section |31 Note, 
though, that in some circumstances extra terms may be needed for the correct analytic 
continuation of the integrals, as discussed in section 4.3.1 of '13'. 

Before concluding this section we would like to mention that while this project was 
in progress and some of the analytic results already obtained as they are written here, 
a preprint by Woynarovich appeared jl5j . Our result 1)5. 24 j) is similar in form to the 
expression proposed by Woynarovich for the 0(1) corrections to the free energy for a 
one dimensional Bose gas with repulsive (5-function interaction. However, there is also 
a major difference. The string of kernels 

+ 02)^02 -03)... - On+l) , (5.31) 

in our (|5.24l) is replaced by a string of the form 

HOi,O2m02, Os)... ^{9n, On+l) , (5.32) 

with 2ip{ki,k2) = (j){ki + k2) + 4>iki - k2) = K{ki,k2) in eq. (3.28) of 15 §. That terms 
of the type (|5.31|) appear in our formulae and not the expression (|5.32|) is unmistakably 
emerging from the Jacobians for the change of variable {rii} {6i} and from their 
definitions as determinants. Woynarovich obtained his result by a calculation of the 
next-to-leading contributions to the free energy, evaluating corrections to the standard 
saddle point result. Such a direct computation would be a highly desirable alternative 
to the more indirect approach taken in this paper. Unfortunately, as stated in the 
paragraph after eq. (5.8) in section V of 15 , for the field theory case the result of |15| 
is divergent in the ultraviolet, R = 1/T 0, limit. This rules out the possibility of its 
consistent agreement with a (^-function defined in (perturbed) conformal field theory, of 
the sort studied in jl() |ll3| [n] and this paper. Nevertheless, the mathematical similarity 
between the final outcomes is striking and deserves further investigation. To make a 
more precise comparison note that 

21og4'^'' P^P^^) ^ [-T-HAF + <Po + cPl) + AS ] (5.33) 

and that —T~^(AF + (j)Q + cpi) matches the first, single-integral, term on the RHS of 
(|5.24|) . AS* should then be compared with the infinite series in (|5.24|) . In eq. (5.8) of 
[15] Woynoarovich notes that his AS can be written as a sum of two contributions: an 
UV convergent part corresponding to 1/2 of our infinite series, plus an UV divergent 
term which has no counterpart in (jniSJ). Thus, in spite of the apparent similarities 
between our results and those of |15j, there are also important discrepancies, which we 
are currently unable to resolve physically. 

^Notice that the kernels in the two cases are actually different, but the derivation in 1151 is quite 
general, and the result is independent of the precise functional form of the kernel. 
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6 Conclusions 



This paper concerned the off-critical version of the boundary entropy g as defined in field 
theory via the identity (|1.6|) . It was shown numericahy that the asymptotic infrared 
expansion for In g, obtained using the Bethe ansatz, matches UV results from conformal 
perturbation theory and the BTCSA at intermediate scales. This was a crucial step in 
the analysis, because it meant that these two alternative definitions are equivalent, and it 
opened up the interesting possibility of deriving an exact expression for the conformally- 
perturbed (7-function by using the Bethe ansatz technique. The first step toward this was 
to give the exact prescription, as the width of the strip R tends to infinity, to transform 
sums over the quantum numbers into integrals on rapidity variables. The idea is that 
the subleading i?-independent terms in this expansion should build up to form the 
boundary entropy g. A careful inspection of this expansion, motivated by the plausible 
assumption that the final result should depend on the bare single-particle energies only 
through their dressed versions, i.e. through the TBA pseudoenergies e{9), led to a partial 
resummation with corrections written purely in terms of £{9)- The expressions for \ng 
written in (|5.24|) and in (|5.30j) are the main new results of the paper. Equation (|5.24j) was 
carefully checked against results obtained using a combination of conformal perturbation 
theory and the Bethe ansatz. The agreement was extremely good, and showed that the 
series is rapidly convergent even in the ultraviolet region (see equations (|5.27p . 1)5. 28() 
and (|5.29|) ). Although it relied at various points on conjectures, we would also like to 
stress that our derivation avoided some of the pitfalls that potentially afflict more direct 
computations of the ^-function: by working in the I — > 00 limit, we always dealt with 
states in which all constituent particles were well-separated, and so the accuracy of the 
Bethe ansatz wavefunctions for high particle density was not an issue. 

There are many open problems related to this project, the first being that the method 
proposed, for all its virtues, is not particularly elegant and a direct approach would be 
desirable for the generalization to more complicated models. It would also be interesting 
to study the corresponding quantities in theories with non-diagonal scattering and in 
systems with massless excitations in the bulk |2fl) . 

Note added: 

There is a numerical error in the third term of the expansion for Ing^ given in eq. 
(5.29) above, which was pointed out to us by Aliosha Zamolodchikov. This is due to 
an inaccuracy in Mathematica's evaluation of the generalised hypergeometric function 
3F2 , which arises in formula (3.14) of ref. [13] for the associated quantity l2- In fact, 
Aliosha Zamolodchikov has found a simplified expression for I2 , as follows: 



(this contrasts with the value —0.083937990... quoted in ref. [13]). Accordingly, the 
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prediction (5.29) above should be corrected to 

In gjil) 0.3214826953191634 + 0.4836924437346968 x A 

- 0.253117570093371858 x + 0.0775 - 0.036 x^ + 0.0195 + . . . (6.2) 

with an even better match to (5.28). We would like to thank Aliosha Zamolodchikov 
for discussions of this point. 
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A TBA and BA in purely elastic scattering models 

In this appendix we summarise the equations relevant to our analysis. 

A.l Periodic boundary conditions: 

The thermodynamic Bethe ansatz equations are jl6j 



€aie) = MaL cosh e - de'cPab{o-e')Lbie'), (a = i,...,iv). (a.i) 



The ground state energy on a circle is expressed in terms of the functions La{9) = 
ln(l + e"^''(^)) as 





(A.2) 



where £Mf L is the bulk contribution to the energy and 



Kb{9) = -^^ln Sab{9). 



(A.3) 
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A. 2 («,/?) boundary conditions: 

The (R-channel) thermodynamic Bethe ansatz equations are jlL)| 

N 

ea{e) = 2MaRcoshe- In fe)(^f " ^)4"H^f + ^)) " / de'(Pab{0 - 0') LbiO') (A.4) 

where a = 1, . . . , N ; the ground state energy on an interval of length R is then 

^ r M 

<"P(M, R) = -y^ -rMa cosh e La{e) + £Ml R + f^ + fp, (A.5) 

where the constant £ is the same as in ()A.2() . and are -R-independent contributions 

to the energy from the boundaries and {R^a\6) , R^p\o)^ are the reflection amplitudes 
corresponding to the two boundary conditions a and /?. Generalisations of these equa- 
tions govern the excited state energies En'^^^iM^R) (12j|^, but in the large R limit we 
are interested in, they reduce to simple (Bethe Ansatz) forms. Suppose that the n^^ 
excited state is made up of m = X^o^i '^^"'^ particles, m^"^ being the number of particles 
of type a. Then 

N m'.'^'l 

'P(M, R) - Ef'^{M, R) = Y,Y.^'' ^i"'^ + 0{e-^^^) , (A.6) 

a=l i=l 

where sums on the RHS with m(") = are understood to be omitted, and the sets of 
numbers {^1"^} satisfy the Bethe ansatz equations 

27rnS") = 2M„/?sinh0;'^)-iln(i?i'^)(0f))i?J^)(e;"))) 

N N 

6=1 jj^i 6=1 j^i 

It is to be noted that the logarithmic branches in the Bethe ansatz equations cause 
some difficulties in the numerics. We impose the branch cut at —it so that the function 
—iln{RR) and any of the terms — iln(— 5) in HA.7|) take values in the range (— 7r,7r]. 
This choice renders the Bethe ansatz ()A.7() fully anti-symmetric (a change of sign in 
any of the quantum numbers n-"'' — > —n^'^^ corresponds to a change ^j"^ — > —O^""^) and 
one can consistently restrict {ra-"^} to strictly positive integers only. 

^Sometimes such generalisations are required even to describe the ground state correctly |12|. but 
these cases will not concern us here. 
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B Some exact results for the boundary Ising model 

The boundary scattering matrix for the free Majorana fermion is [H] 

/ ivr 9\ k — i sinh 9 
Rk(e) =itanh{ — - -] . . , B.l 

where k = 1 — and h is the boundary magnetic field. Therefore, we have 

1/1 4A;cosh6' \ 

"^'^^^ - n [-^e " cosh(20) + 2P - 1 ) ■ ^^-2) 

In this appendix we would like to report exact expressions for 

In^free = ln<7fc=i(/) = 'l jj'^ (^^ + ^^^) + e"' (B.3) 

and 

In^flxed = ln5fe=-oo(0 = -\j/^ ('^(^) - ^^^^) ^^^(1 + • (B.4) 

These are obtained from the following identity, which can be deduced by studying the 
monodromies of the integral (|2.8|) along the lines sketched in J21': 

/ de \- In (1 + e-"™^^^'') = ln(2) - -(1 + ln(^/x) - 7i?) - 2S{x) (B.5) 

vr cosh vr 

where 7E = 0.57721566 ... is the Euler-Mascheroni constant and 

oo / 

S{x) = J] In 



n=l 



+ (2n - l)vr - Y^(2n - l)2vr2 + 



(2n - l)vr + v^(2n - l)2vr2 + x2 



The idea of the derivation is to determine the positions of the singularities in ()B.6|) 
using the pinched-contour argument of [2^, adding counterterms to make the infinite 
sum S{x) convergent. The remaining parts of HB.5|) were then fixed by studying the 
X ^ limit. (Alternatively, (|B.5|) can be proved using the Bessel- function technique of 
[mi^ and an identity due to Schlomilch [211 •) From HB.5|1 . we have 

41n5free/flxed(0 = " ln(l + e"') T (ln(2) - ^(1 + ln(^//) - 7^5) - 2S{1)^ . (B.7) 
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